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ABSTRACT 



We develop a new approach to the well-studied anti-correlation between the 
optical-to- X-ray spectral index, aox, and the monochromatic optical luminosity, 
lopt- By cross-correlating the SDSS DR5 quasar catalog with the XMM-Newton 
archive, we create a sample of 327 quasars with X-ray S/N > 6, where both 
optical and X-ray spectra are available. This allows aox to be defined at arbitrary 
frequencies, rather than the standard 2500 A and 2 keV. We find that while the 
choice of optical wavelength does not strongly infiuence the aox — hpt relation, the 
slope of the relation does depend on the choice of X-ray energy. The slope of the 
relation becomes steeper when aox is defined at low (~ 1 keV) X-ray energies. 
This change is significant when compared to the slope predicted by a decrease in 
the baseline over which aox is defined. The slopes are also marginally fiatter than 
predicted at high (~ 10 keV) X-ray energies. Partial correlation tests show that 
while the primary driver of aox is lopt, the Eddington ratio correlates strongly 
with aox when lopt is taken into account, so accretion rate may help explain these 
results. We combine the aox — lopt and T-Lhoi/l-'Edd relations to natu rally explain 



two re sults: 1) the existence of the T — Ix relation as reported in lYoung et al. 



(120091 ) and 2) the lack of a F — lopt relation. The consistency of the optical/X-ray 
correlations establishes a more complete framework for understanding the relation 
between quasar emission mechanisms. We also discuss two correlations with the 
hard X-ray bolometric correction, which we show correlates with both aox and 
Eddington ratio. This confirms that an increase in accretion rate correlates with 
a decrease in the fraction of up-scattered disk photons. 
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Introduction 



Optical/UV and X-ray continua in quasars are tliougtit to originate in two physically dis- 
tinct, but nevertheless related r egions. The optica l /UV spectrum is dominated by a Big Blue 



Bump continuum featur e fe.g.. lElvis et al 



from an accretion disk (jShields et al 



1978 



19941). most likely from thermal emission ar ising 



Malkan fc Sargentlll982l : IWard et all 119871 1 A 



hot gas lying above the disk in some unknow n geometry then up-scatters these disk p hotons 



to form a hard p ower-law (< >~ -1) ('e.g..lMateos et al. 



X-ray spectrum (IHaardt fc Maraschilll991 



Zdziarski et al. 



2005 



2000 



Mainieri et al. 


20 


07 


) in 


Kawaguchi et al. 


2001) 



While the shapes of the optical and X-ray continua are well known, the connecting 
continuum in the extreme ultraviolet (EUV) is not, because of absorption from our own 
galaxy. The EUV and soft X-rays form the ionizing continuum (hz/ > 13.6 eV) needed to 
ionize a quasar's characteristic broad lines. The ionizing cont inuum is also responsible for 
any line-driven outflows from the accretion disk (jProgal 120071 ). Since direct observation of 
this continuum is not possible, it is often parameterized by an imaginary power-law with 
spectral index running from 2500 A to 2 keV in the rest frame. 

The importance of the ionizing c ontinuum motivates the decades long study of the 



anti-correlation between and Lq^oo (lAvni fc Tananbaum 



Tananbaum et al 



1982 



19861: lAnderson fc Mareonl 1 19871 : 1 Wilkes et al.lll994l : IPickering et al.iil994 



Kriss fc Canazares 1985: 



Avni. Worrall fc Margon Il995[ ): however, controversy continues. While so me studies have 



found that ana-, depends primarily on redshift rather than optical luminosity (IBechtold et al. 



20031 : iKelly et al.l 120071 ) , others have suggested that an^. may be entirely independent of 



redshift and optical luminosity (lYuan et al.lll998l : iTang et al.l 120071 ). arguing that if the dis- 
persion in L2500 is much greater than the dispersion in \jikev^ an apparent but artificial 
correlation will result. A narrow range in optical luminosity will enhance this effect. 

In recent years, studies of the — -^2500 relation have become more comprehensive, 
with higher X-ray detection rates and largely complete samples that span a wide range o f 



luminosities ( Vignali et al. 
For example. 



Steffen et al. 



2OO3I: Istrateva et al.l[2005l : Isteffen et allbood : I Just et al.l[2007[ ). 



(I2OO6I ) work spans five and four decades in UV and X-ray lumi- 
nosity, respectively, which helps mii iimize the effec t of a n artificial aox — -^2500 correlation 
due to selection effects. In addition, IStrateva et al.l (120051 ) estimated the dispersion in L2500 
and \j2kev-, finding that the difference in dispersions was not enough to artificially reproduce 
a correlation. These studies used partial correlation analysis to conclude that the aox— z 
correlation is insignificant, while the aox — -^^2500 anti-correlation is highly significant. How- 
ever, while the significance of the relation continues to b e better understood, little is known 



of the underlying phy s ics. Current disk-corona models (iBelo 



1999: 


Navakshin 


2000: 


Malzac 


2001: 



jorodov 1999: Zdziarski et al. 



b|) do not directly predict the 
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Q-ox - L2500 relation. 

Previous studies of the aox — -^^2500 relation have largely lacked spectral slopes in the 
X-rays, and so used the traditional endpoints of 2500 A and 2 keV, where an X-ray spectral 
slope (r ~ 2, where T = -a + 1 for Fj, cx i/") is assumed in order to obtain the X-ray flux at 
2 keV. A systematic study with both optical and X-ray spectra enables an investigation of 
the relation at different frequencies than those traditionally used, revealing clues about the 
relation's physical underpinnings. 

The sample selection and measurement of aox are described in §2, and results of the data 
analysis are given in §3. Section 4 discusses the results, specifically their connection to other 
optical/X-ray relations found for quasars. We assume a standar d cosmology throug hout the 
paper, where Hq = 70 km s~^ Mpc~^, = 0.3, and = 0.7 (jSpergel et al. 2003). 



2. The Sample 



We have cross-correlated the DR5 Sloan Digital Sky Survey (SDSS) with XMM-Newton 
archival observations to obtain 792 X-ray observations of SDSS quasars. This gives the largest 
sample of opti cally selected quasar s with both optical and X-ray spectra (473 quasars with X- 
ray S/N > 6). lYoung et al.l (120091 ) describes the SDSS quasar selection and cross-correlation 
with the XMM-Newton archive, as well as the X-ray data reduction and spectral fitting. For 
the purposes of this paper, we exclude broad absorption line (BAL) and radio-loud (RL) 
quasars, as th ese are known to have different a„„. distr ibutions than typical, radio-quiet 



(RQ) quasars (iGreen fc Mathui 



an incomplete BAL list (IShen et al. 



the DR5 BAL catalog (I Gibson et al. 



19961: lYoung et al.ll2009f ). While lYoung et all fl2009f ) used 
20081) , this paper updates the BAL classification using 



20091 ). A total of 18 sources are re-classified as BALs, 



and an additional 15 sources are re-classified as non-BALs. A net result of three additional 
BALs, gives the SDSS/XMM-Newton Quasar Survey 55 BALs in total. 

Since the goal is to use our knowledge of the X-ray spectrum to determine the effects 
of the X-ray continuum shape on the calculation of aox, we form a sample of quasars with 
enough X-ray S/N to fit a spectr um (S/N > 6 ) in th e primary sample of this study, hereafter 
referred to as SPECTRA. As in lYoung et al.l (120091 ). we fit three models to each spectrum: 
1) a single power-law (SPL) with no intrinsic absorption, 2) a fixed power-law (FPL) with 
intrinsic absorption left free to vary, and 3) an intrinsically absorbed power- law (APL), with 
both photon index P and absorption Nh left free to vary. Any spectrum that does not give 
a good fit with any of the above models is excluded from the sample, in order to ensure 
well-determined X-ray fluxes. Therefore the final sample will not include any sources with. 
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for example, a significant contribution from a strong soft excess component. In addition, 
any spectrum tliat prefers eitfier of the absorbed models (FPL or APL) is excluded from 
the sample, in order to minimize the effect of absorption on both X-ray and optical flux 
measurements. This reduces the sample by 55 sources, or 11.6% of quasars with X-ray S/N 
> 6. Once all of the selections are taken into account, the sample contains 327 quasars 
spanning a redshift range of z = 0.1 - 4.4. 

The median upper-limit for intrinsic absorption in the SPECTRA sample is Nj:/ < 4 x 
10^^ cm~^. To check for the effect of undetected absorption in low S/N X-ray spectra, we 
raise the threshold to S/N > 20, where Nh upper limits decrease to ~ 6 x lO^'' cm~^. This 
sample contains 99 sources and will hereafter be referred to as HIGHSNR. 

Finally, we avoid bias against X-ray weak sources by adding undetected or low-signal 
sources to the SPECTRA sample, to create the CENSORED sample. To include sources 
with S/N < 6, we apply a fixed power-law fit with no absorption, where the power-law is 



fixed to r = 1.9, the average value for the SDSS/XMM- Newton Quasar Survey (lYoung et al. 



20091 ) . RL and BAL quasars are still excluded, as are any sources with bad X-ray fits, leaving 



514 sources in the sample spanning a redshift range of z = 0.1 - 5.0. 

The three samples used in this paper are described in Table 1. Figure [1] gives the i band 
magnitude vs. redshift distribution of each of the three samples. 



2.1. Measuring Oox 

For each SPECTRA sample source, we fit six versions of an unabsorbed power-law. 
The normalization is frozen to rest-frame energy values of 1, 1.5, 2, 4, 7 and 10 keV and 
then the monochromatic fluxes are taken at each energy. Freezing the normalization at each 
energy enables us to accurately calculate the error on the flux based on the error on the 
normalization of the fit. If the normalization was frozen to 1 keV observed-frame, as is 
standard, both the error in normalization and the error on F would affect the error on the 
flux. For 215 quasars, at least one of the rest-frame energies lies below the observed range. 
In these cases, we normalize the spectrum to 1 keV and extrapolate to obtain the flux at the 
rest-frame energy. 

We do not normalize to rest-frame energies lower than 1 keV because undetected ab- 
sorption with a median upper limit of Nh ~ 5 x 10^^ would eliminate any measured flux at 
0.5 keV, and would reduce the flux measured at 0.75 keV by an order of magnitude. At 1 
keV, the change in flux is only a factor of ~2, which would change by ~0.1. Limiting 
the normalization energies to 1 keV and higher has the added benefit of minimizing the 
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contribution from any soft excess component. 

We also do not attempt to fit fluxes at rest-frame energies higher than 10 keV, as 
these are only directly available for high redshift (and therefore high luminosity) sources. 
Extrapolation with the fitted power-law for low redshift sources would not take into account 
the reflection bump that begins rising at 10 keV and peaks near 30-50 keV. For low redshift 
sources, rest-frame energies of 7-10 keV will sample the 7-10 keV observed spectrum, where 
under- or over-subtraction of the background can alter the measured fluxes. The resulting 
increase in flux measurement error will affect the dispersion of the — ^opt relation, which 
will be discussed further in §3. 

Before calculating optical luminosities, we first correct the magnitudes for Galactic red- 
dening. Corrections are applied using a Galactic dust-to-gas ratio and the Galactic hydrogen 
column density from WebPIMM^, which is based on the 21 cm HI compilation of Dickey & 
Lockman (1990) and Kalberla et al. (2005). We then interpolate a power-law through the 
well-calibrated SDSS photometry (~2-3% photometric error). As detailed below, we care- 
fully select wavebands for interpolation to avoid host galaxy contribution at low redshift, 
the Lya forest at high redshift, and major emission lines. 

For redshifts z < 0.55, the host galaxy can redden the magnitudes at long wavelengths, 
so we use the bluest bands, u, g and r, to determine the optical slope, while also avoiding 
those bands containing the Mgll line. For z < 0.25, we use the g and r, switching to u and 
g for 0.25 < z < 0.43, and then to u and r for 0.43 < z < 0.55. (Note that the Mgll line lies 
largely between the u and g filters for the middle redshift range, 0.25 < z < 0.43.) 

For quasars at redshifts above z = 0.55, we assume host galaxy contribution to be 



minimal (IRichards et al.ll2003l ) and fit a power-law through all wavebands using an ordinary 
least-squares regression. While lines such as Mgll, GUI and GIV affect the photometry at 
various redshifts, their addition to the continuum flux is minimized by the broad bandwidths 
(the g,r,i and z bands have FWHM ~ 1000 A). In addition, the regression fit to all wavebands 
minimizes the contribution of any one waveband that is an outlier due to line emission. 

As the Lya line begins to affect the wavebands, they are removed from the regression 
fit. The u band is removed from the fit at z = 1.8, and the g band at z = 2.3. At z > 3.75, 
for two quasars, the Lya line af fects the spectrum s i gnific antly and we simply fix the optical 



slope to the median value from IVanden Berk et al.l (120011 ) . aopt = 0.46. 



The normalization of the optical power-law is also adjusted according to redshift. For 
sources with z < 0.55, the power-law is normalized at the bluest band to avoid host galaxy 



^http: / /heasarc. gsfc.nasa.gov/Tools/w3pinims. html 
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contamination. For sources with 0.55 < z < 3.75, the power-law is normahzed using the 
ordinary least squares fit to all of the available magnitudes. For the single source with z > 
3.75, the power-law is normalized to the z band to minimize the effect of the Lya forest. 
From the measurement of the index and normalization of the optical power-law, we can 
extrapolate to three monochromatic optical luminosities at 1500, 2500 and 5000 A, where 
the log of the optical luminosity is hereafter denoted as lo^t- 

The blend of Fell and Balmer emission lines that make up the S mall Blue Bump 



(IWills. Netzer fc Wills Ill985l ) affects monochromatic optical fluxes by 0-45% (IBoroson fc Green 
19921 ). This causes at most a small overestimate of aox by A(aoa;) ~ 0.06. 



The mean optical slope obta ined for the SPSS /XMM- Newton Quasar Survey [o-ovt = 
-0.40) is similar to the slope of the lVanden Berk et al.l (120011 ) mean composite (aop^ = -0.44), 
indicating that the host galaxy and emission line contributions have been adequately taken 
into account. The di spersion of the optica l slope is Oa = 0.30, similar to the dispersion seen 
in the composites of iRichards et al.l (120031 ) . where cTq, ~ 0.26. 



Once the optical luminosities are obtained, corrections for dust-reddening are applied 
where necessary. We obtain an estimate of E(B-V) from the relative {g — i) color, where 
relative colors compare a quasar's measure d colors with the med ian colors in its redshift bin, 
so that A((7 — i) = {g — i) - <{g — i)>z (IRichards et al.ll2003l ). Quasars with A{g — i) < 
0.2 lie within the normal color distribution for SDSS quasars, and so are not likely to be 
affected by dust-reddening. For quasars with A{g — i) > 0.2 the E(B-V) estimated from 
the rel ative colors is used to de-redden the optical luminosities with an SMC-type extinction 
curve (jPrevot et al.l Il984j ) . Note that the correction can only be made for sources with 
redshift less than 2.2 for the relative colors to be valid, due to the entry of the Lya line 
into the spectrum. Sources with E(B-V) > 0.04 are classified as dust-reddened and make up 
7.6% of the SDSS/XM M-Newton Quasar Su rvey. This fraction is comparable to that found 
for SDSS quasars 



Richards et al.ll2003l ). 



We next calculate aox with the standard formula: 

log{fx/fo) 



an 



log{i^x/i^o) 

where is the monochromatic flux at 1, 1.5, 2, 4, 7, or 10 keV and fo is the monochromatic 
flux at 1500, 2500, or 5000 A. 

For each aox — lopt relation, we fit an ordinary least-squares regression line and obtain 
the slope and intercept. The errors are at the la level, and the standard deviation on the 
regression gives the observed scatter in the relation. 



For the CENSORED sample, we fix the power-law to the sample average of F = 1.9 for 
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187 sources with S/N < 6, and then fit the monochromatic X-ray flux. For 59 undetected 
sources with S/N < 2, we flt the 90% upper-hmit on the flux. Since the data now include 
upper hmits, we input the hsts of and lo^t into the 'Estimate and Maximize' (EM) 
function to compute the hnear regression. This m ethod is among sever al routines available 
on ASURV Rev 1.2, a su rvival statist i cs pa ckage (jLavalley et al.l 119921 ). which implements 
the methods presented in llsobe et al.l (119861 ). Survival statistics account for censored data 
by assuming that the upper limits follow the same distribution as the detected points. The 
EM regression reduces to an ordinary least squares flt for uncensored data. 

For all three samples, we use a partial correlation test designed to measure the corre- 
lation between two variables, while controlling for the effects of a third. The FORTRAN 
program CENS_TAU calculates Kendall's ri2,3 partial correlation coefficient, and can tak e 
censored data into account using the methodology presented in (lAkritas fc SiebertI Il996l ). 
Unless otherwise noted, redshift is the third variable for the correlation tests. Both ASURV 
and CENS.TAU are made available by the Penn State Center tot ABtto^tatiBticfl 



Results 



We first measure the ft nrr — Irmt anti - correlation using th e analysis of previous studies 



dVignah et aP I2OO3I: Istrateva et aLllpOOsI : Isteffen et all bood ) on the CENSORED sample. 
In comparison to Steffen et al. (j2006l ). the CENSORED sample is larger (N=514, compared 
to N=333), with a similar range of ~4 decades in X-ra y luminosity. The X -ray upper limit 
fraction is similar as well, (11%, compa red to 12% for ISteffen et al.l (120061 )). The range in 
optical luminosity in lSteffen et al.l (120061 ) is roughly 1.5 decades larger than the range in this 
paper (3.5 decades), due largely to their inclusion of low and moderate-luminosity samples 
in addition to the SDSS. 

Figure [2] shows the — Lpt anti-correlation, significant at the 9.0cr level (ti2,3 = -0.20) 
in the SDSS/XMM Quasar Survey when defined at the conventional 2500 A and 2 keV 
fiducial points. The solid line is the best-fit EM regression line for the CENSORED sample: 

aox = (2.133 ± 0.388) - (0.122 ± 0.013)/o^(L ») 



The slope is consistent within errors to the best-fi t regression from ISteffen et al.l (120061 ) 
(dashed line in Fig. [2]). Unlike in lSteffen et al.l (120061 ). we are not able to reject a correlation 



between 



and redshift while controlling for L ° . Instead we find the partial correlation 

° 2500A 



^http: / /www. astrostatistics.psu.edu/statcodes/cens_tau 



- 8 - 



to be marginally significant at the 3.8cr level (ti2,3 = -0.085) in the CENSORED sample. 
This is likely due to the smaller range in optical luminosity, as well as the large scatter in 
aox- Note that in the SPECTRA sample, the partial correlation between and redshift is 
largely insignificant, except at 1 keV where, again, the scatter is significant (see Tabled]). 

We can now test the correlations of aox versus optical luminosity using the SPECTRA 
sample, for each of the six X-ray energies and three optical wavelengths where Uox is defined. 
Figure 3 shows 6 examples of the relation, defined at 1500, 2500 and 5000 A and at 1 and 
10 keV. The plots show a trend in the slope and dispersion of ttox — ^opt relation: the slope 
fiattens and the dispersion tightens as the X-ray energy used to define aox increases. It is 
clear from the plots that the effect of changing X-ray energy is much stronger than that of 
changing optical wavelength. 

The definition of aox depends on the baseline over which it is defined; lengthening the 
baseline can produce an artificial fiattening of the aox — hpt slope. To check that the intrinsic 
slope of the correlation is changing, we define the standard slope by the relation at 2500 A 
and 2 keV. We then calculate how the change in baseline affects the slope and intercept of 
the relation. As the changing baseline frequencies are accompanied by a change in the X-ray 
flux, we include this effect by assuming an average value for the X-ray slope (< ax >= —0.9). 
The equations for the change in slope (m) and intercept (b) then become: 

m = m ^^^^^ 



and 

h' 



logiPx/vo) 



logiux/iyo) 

where Ux is the X-ray frequency at 1, 1.5, 2, 4, 7, or 10 keV and Uo is the optical frequency 
at 1500, 2500, or 5000 A. 

A change in the optical baseline, for example from 2500 to 5000 A, affects the slope 
of the anx — Inpt relati on by an amount comparable to the error. This is in contrast to 



Vasudevan et al.l (120091 ) . who flnd a marginal strengthening of the aox — Lpt relation as the 
optical reference point is moved to shorter wavelengths. This may be due to a steepening 
of the correlation slope, but the small sample size and the strong effect of reddening on the 
sample do not allow a deflnitive statement. Since the slope does not depend signiflcantly on 
optical wavelength in our results, we plot our results only for 2500 A in Figure 4, in order 
to focus on the change with X-ray energy. 



Figure 4a shows that while changing the X-ray baseline clearly has an influence on the 
measured slopes, the slope for 1 keV is signiflcantly steeper than it would be if the intrinsic 
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correlation were constant. The slopes at 4, 7 and 10 keV are flatter than predicted by the 
baseline effect at the la level. However, since the baseline effect is arbitrarily normalized 
to the slope at 2 keV, the main result of subtracting the baseline effect is not to assess the 
significance of the slope at a given energy, but rather to show that the intrinsic slope of the 
c^ox — lopt relation is flattening as X-ray energy increases. 

Figure 46 and 4c show the intercept and dispersion of each relation as measured via 
the least squares flt. The intercept of the relation decreases signiflcantly with X-ray energy 
compared to the predicted curve. The dispersion of the — hpt relation largely follows the 
predicted curve, though it is larger than predicted at the lowest and highest energies. This 
is likely due to the larger errors in X-ray flux at the ends of the spectrum. 

For low redshift sources, rest-frame energies of 7-10 keV will sample the 7-10 keV ob- 
served spectrum, where under- or over-subtraction of the background can alter the measured 
fluxes and increase the error in the flux measurement. When we exclude sources with z < 
0.5, we flnd that the dispersion of the aox — lopt relation at 7 and 10 keV decreases by ~ 
0.005, though the dispersions do not decrease all the way to the baseline level. This suggests 
that the background does have some effect at low redshifts, but it is not the only source of 
scatter at high energies. 

Interestingly, the trend in slope and dispersion of the aox — lopt relation for the CEN- 
SORED sample is mostly consistent with the baseline effect (Fig. [S]), with a signiflcantly 
steeper slope than expected at 1 keV. The slopes also show the same marginal flattening 
at high X-ray energies as the SPECTRA sample. These results are not as strong as in the 
SPECTRA sample. 

Figure [2] suggests that the main effect of including low S/N and undetected sources is an 
increase in scatter, diluting any potential trends. The red circles in the flgure, representing 
detected sources with no spectral flt, show a weaker anti-correlation, while the blue circles, 
representing sources with spectral flts, are more strongly anti-correlated. A partial correla- 
tion test shows that the correlation decreases in signiflcance from 8.6cr (ti2,3 = -0.284) for 
sources with S/N > 6 to 4.6cr (ri2,3 = -0.131) for sources with S/N < 6. 

A higher percentage of S/N < 6 sources are X-ray weak {aox < -1-8): 18.7%, compared to 
only 3.7% in the SPECTRA (S/N > 6) sample, where absorbed sources have been excluded. 
Since absorption cannot be accounted for in the low S/N sources, it is a likely source of the 
scatter. Figure E] conflrms that though the SPECTRA sample does not remove every X-ray 
weak source, those with the heaviest optical obscuration {/S.{g — «) > 0.3) are excluded from 
the sample. This illustrates the importance of X-ray spectra in excluding absorbed sources. 

For the HIGHSNR sample, the slopes show the same trend as found in the SPECTRA 
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sample (cf. Figures 4a and 6a), though the errors are larger due to the smaller sample size. 
Therefore the significant change in slope of the — hpt relation is not likely to be due to 
the effect of undetected absorption in the sources with low S/N X-ray spectra. However, 
the scatter in the relation defined at low energies does decrease for the HIGHSNR sample 
(cf. Figures 4c and 6b), suggesting that undetected absorption affects the dispersion of the 
Oiox — lopt relation when sources with low S/N spectra are included. 

Results for the CENSORED, SPECTRA and HIGHSNR samples are summarized in 
Table [H which gives the slope, intercept and dispersion for the — lopt relation at each 
X-ray energy for each sample. 



4. Discussion 



The SPECTRA sample will be used for the following discussion. 



4.1. Partial Correlations with Physical Parameters 



To gain increased understanding of the trends found above, we tested for correlations 
with physical parameters. Controlling for the effects of the optical luminosity, we computed 
the partial correlation coefficient between aox and black hole mass, Eddington ratio, and 
redshift. 



Black hole masses were derived in IShen et al.l (120081 ) using the width and continuum 
luminosities of broad emission lines in the optical spectrum: H/3 for < z < 0.7, Mgll for 0.7 
< 1.9, and CIV for 1.9 < z < 4.5. The broad line region (BLR) is assumed to be virialized, 
so that the continuum luminosities are related to the BLR radius by the R-L scaling relation, 
and the emission line widths give the velocities of the BLR clouds. Further details o f the 
line-fitting and scaling relations used for e ach line can be found in IShen et al.l (120081 ). To 
calculate the Eddington ratio, we use the [Richards et al.l (120061 ) bolometric correction for 
3000 A, K.snnn = 5.62i:1.14 . Though bolometric corrections are by nature uncertain. Figure 
12 of [Richards et al.l (120061 ) shows that 3000, as well as 5100 A are regions of relatively small 
dispersion in the composite SEDs. 



Tables [3ll5] list the partial correlation coefficients (ri2.3) and their significance level at 
each frequency for which aox is defined. Table |3] compares aox with black hole mass (M^//), 
Table m compares aox with the Eddington ratio {Lboi/^Edd)^ and Table O compares aox with 
redshift (z). We find that the primary dependence is of aox on optical luminosity. However, 
there is a significant partial correlation with Eddington ratio that increases with X-ray 
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energy, with no dependence on optical wavelength (see Fig [8]). The strength of the partial 
correlation suggests that Eddington ratio may play a role in determining aox- 

The partial correlation tests show a marginally significant partial correlation between 
and redshift (Table [5]) but the correlation disappears at higher energies. This is likely 
due to the large scatter in the aox — Lopt relation at 1 keV. The tests show no significant 
partial correlation between aox and Mbh (Table [3]). The black hole masses of the SPECTRA 
sample range from log Mbh = 7.4 to 10 Mq. 

Note that Kelly et al. (2007; see their Fig. 8) found problems with the CENS_TAU 
program, finding that zero partial correlation between two variables does not necessarily 
result in a null correlation coefficient (r = 0). Assessing the real significance for a given 
value of the coefficient is therefore difficult at low values. However, given the high value of 
the partial correlation between aox and Eddington ratio at high X-ray energies (ti2,3 ~ 0.3), 
a significant correlation appears likely. 



4.2. Consistency Between aox — ^opti F — Ix and Y—Vii,oil'^Edd Relations 
In this s ection , we compare the aox ~ ^opt relation to the T — relation found in 



Young et al.l (120091 ). An apparent contradiction arises, as these relations should lead to 



a secondary relation between F and lopt, when in fact the F — lopt relation is found to be 
insignificant in the SDSS/XMM- Newton Quasar Survey. This contradiction can be resolved 
by accounting for the dispersion in the primary relations. 



The relation between F and l^ (lYoung et al.ll2009l ) depends on the X-ray energy. The 
correlation is marginally significant and positive at the lowest energy {x = 0.7 keV) and 
highly significant and negative at higher energies [x = 4-20 keV). The change in slope of 
the correlation indicates a pivot in the X-ray spectrum near 1 keV, where the correlation is 
insignificant. 

At 10 keV, the T — l^ relation is described by 

r = -(0.335 ± 0.032)/iofceV/ + (10.732 ± 0.818), a = 0.31, (1) 
The aox — Lpt correlation at 2500 A and 10 keV is described by 

= -(0.081 ± 0.010)/25oo + (1-044 ± 0.319), a = 0.11, (2) 

which is equivalent to 

Ix = (0.734 ± 0.034)/25oo + (3.449 ± 1.054), a = 0.37. (3) 
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Combining equations ([T]) and ([3]) should give the following correlation between F and lopt'- 

r = -(0.246 ± 0.0 26)/25oo + (11.887 ± 0.898) 

Such a correlation is not obseryed to be significant in the SDSS/XMM-Newton Survey 
(lYoung et al.ll2009l : iRisaliti et al.l 120091 ). so the relations appear to be inconsistent. 



This apparent contradiction can be explained by accounting for the dispersion associ- 
ated with the T—Lboi/LEdd relation. The hard (2-10 k eV) X-ray slope ha s been shown to 



2006; 


Kellv et al. 


2007; 


Risaliti et al. 


r vv J 

2009). 



by analyzing a sample of 35 quasars over three orders of magnitude in luminosity, finding 
a stronger T-Uni/L EdA than r-FWHM(H/3). A linear fit to the SDSS/XMM-Newton data 
fiRisahti et aLll2009h gives: 



r = (0.31 ± OM)log{Lhoi/LEdd) + (1-97 ± 0.02), a = 0.33 



(4) 



To illustrate the discussion below, a notional SED is shown in Figure M The figure 
shows a simple model of a quasar SED, where specific luminosities are measured at 2500 
A and 1 keV. Since the T — relation is not significant at 1 keV, using the luminosity 
at this energy (instead of at 10 keV, as above) allows the focus to be on the effect of the 
dispersion associated with the T—Lhoi/LEdd relation. In the figure, model (A) illustrates a 
source with lower lopt (10^^ ergs s~^) and model (B) shows a source with higher lopt (10^^ 
ergs s"^). Inserting the optical luminosity into the appropriate aox — lopt relations gives an 
average value for aox (dashed lines) and therefore for hkev and hokev- 



Since 1 keV is near the pivot point of the X-ray spectrum (lYoung et al.l 120091 ). the 
average X-ray slope will not depend on likeVy and so will just be the mean F of the sample 
(< F >~ 1.9). The dispersion around the mean will arise from two sources: 1) dispersion 
associated with the T—Lf,oi/LEdd relation (Eqn HI shown in Figure [9]) and 2) dispersion in 
Lboi/ljEdd for a given lopt- While the fiattening of the aox — lopt relation with energy implies 
harder X-ray slope as optical luminosity increases. Fig. O shows that this expected hardening 
is small compared to the dispersion around the mean. Therefore, any secondary relation 
between F — lopt would have too large a dispersion for the correlation to be significant. 

Moreover, Figure [9] shows how the F — relation results naturally as a secondary 
correlation with large scatter. The figure shows that sources at a given lopt with higher 
hard X-ray luminosities will have harder X-ray slopes, while those with higher soft X-ray 
luminosities will have softer slopes. The combination of the aox — lopt relation with the F- 
Lboi/l-'Edd relation thus returns the T — Ix relations for X-ray energies a; < 1 keV or x > 1 
keV. 
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4.3. Bolometric Corrections 

A s Qio.r can be used t o apply a bolometric correction to the hard X-rays (2-10 keV) 



[e.g., 



Hopkins et al.l 120071 ) by providing a normalization for the X-ray spectrum, we plot 
in Figure [To] the dependence of K2-iokeV on aox- Since L^oi is tightly related to the optical 
luminosity through the bolometric correction, and L2keV is tightly related to L2_ioA:ey through 
integration over F, the linear correlation is expected. However, both the correlation and its 
dispersion, a = 0.12, are of interest when determining X-ray bolometric corrections. For 
reference, the best-fit line is: 

H2-wkeV = -(2.282 ± 0M5)aox - (1-67 ± 0.10) 

The X-ray weak quasars with very large bolometric corrections (^2-10^6^ > 1000) do not 
have enough S/N for X-ray spectral fits. One object is optically red {A{g — i) = 0.33), so it 
is possible that these extreme objects are suffering from absorption, affecting both their aox 
and K2-iokeV values. An additional factor to consider is accretion rate. We find a significant 
partial correlation between /t2-iofceV and Lboi/ljEdd- The presence of Lboi on both axes can 
result in a false correlation, but a partial correlation test with L^o/ as the third variable shows 
that the correlation is significant at the 7.7a level, with a ri2,3 = 0.232. The correlation is 
plotted in Figure [TT] with a best fit line: 

logK2-iokeV = (0.443 ± 0.058)logLtoi/ LEdd + (2.227 ± 0.045) 



This correlation is consistent with previous results (e.g., IVasudevan fc Fabian! l2009l ). and 
suggests that the fraction of disk photons up-scattered to the X-rays decreases as accretion 
rate increases. 



5. Conclusions 

We find that the slope of the aox — lopt relation depends on the X-ray energy at which 
aox is defined, though there is no significant dependence on optical wavelength. While 
a change in slope is predicted when the baseline over which aox is defined increases or 
decreases, the slopes at the lowest X-ray energy are significantly steeper than predicted, 
while the slopes at high X-ray energies are marginally flatter than predicted. This suggests 
that the efficiency of low energy X-ray photon production depends more strongly on the 
optical/UV photon supply, whereas the efficiency of high energy X-ray photon production 
remains relatively constant with respect to the seed photon supply. Partial correlation tests 
show that Eddington ratio is highly correlated with aox when controlling for lopt^, and this 
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partial correlation increases in significance with X-ray energy, suggesting that accretion rate 
hkely plays a role in the results. Nevertheless, lopt remains the primary driver of aox- 

To check the results, we redo the calculations for a sample including low S/N sources 
(CENSORED) and for a sample including only high S/N sources (HIGHSNR). The CEN- 
SORED sample shows the same trends as in the SPECTRA sample, though the results are 
less significant. This is likely due to the diluting effect of absorption in the low S/N sources, 
which emphasizes the importance of X-ray spectra in eliminating sources with absorption. 
The HIGHSNR sample reproduces the trend in the slope of the aox — lopt relation found for 
the SPECTRA sample, indicating that the results are not affected by undetected absorp- 
tion. The dispersion of the relation does decrease significantly for the HIGHSNR sample, 
suggesting that undetected absorption does introduce some error into the determination of 
X-ray fluxes for sources with moderate S/N. 

By combining the aox — lopt relation with the V-hboi/ljEdd relatio n, we can na t urally 



explain two results: 1) the existence of the F — Ix relation as reported in lYoung et al.l (120091 ) 
and 2) the lack of a F — lopt relation. Given the measured dispersions and correlation 
coefficients, it is likely that the aox — lopt and F- l^hoi/^EM relations are the primary relations, 
while the T — Ix relation is secondary. Understanding the relations in this order will give a 
more complete framework for understanding the relation between optical and X-ray emission 
in quasars. In order to fully understand the physical implications of the results shown here, 
a study of the selection, orientation and variability effects is needed. Such a study is now in 
progress. 

Finally, for the purposes of aiding studies of the hard X-ray bolometric correction, we 
include a description of the 1^2-iokeV - «oz correlation and dispersion. We find a significant 
partial correlation between 1^2-wkeV and hi^oi/^Eddi when taking bolometric luminosity into 
account as the third variable. This is consistent with a decreasing fraction of up-scattered 
X-ray photons as accretion rate increases. 
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Table 1. Sample Definitions 



Sample name l^tot S/N median N^/ z i mag % targets 



CENSORED 514 all 0.39 0.1-5.0 15.2-20.7 -2% 

SPECTRA 327 S/N > 6 0.39 0.1-4.4 15.2-20.4 -3% 

HIGHSNR 99 S/N > 20 0.06 0.1-3.0 15.6-19.9 ~9% 



Note. — All three samples exclude RL and BAL quasars, sources with significant 
intrinsic absorption, and sources with bad fits to the X-ray spectrum. 



Table 2. Results for the — loptRdation 



X-ray Energy 




CENSORED 






SPECTRA 






HIGHSNR 






Slope 


lut orcc^])! 


n 


Slopc^ 


lut orcc^])! 


a 




Iiitcrct^pt 


a 


1 keV 


-0.17±0.02 


3.47±0.47 


0.19 


-0.19±0.02 


4.13±0.46 


0.15 


-0.13±0.02 


2.45±0.54 


0.11 


1.5 keV 


-0.13±0.01 


2.43±0.42 


0.17 


-0.14i0.01 


2.70±0.38 


0.13 


-0.10±0.02 


1.55±0.50 


0.10 


2 keV 


-0.12±0.01 


2.13±0.39 


0.16 


-0.12±0.01 


2.27±0.35 


0.11 


-0.08±.02 


1.09±0.49 


0.10 


4 keV 


-0.099±0.01 


1.49±0.35 


0.14 


-0.098±0.01 


1.54±0.30 


0.10 


-0.06±0.02 


0.30±0.53 


0.11 


7keV 


-0.090±0.01 


1.26±0.34 


0.14 


-0.086±0.01 


1.21±0.31 


0.10 


-0.04±0.02 


-0.15±0.58 


0.12 


10 keV 


-0.080±0.01 


0.97±0.34 


0.13 


-0.081ib0.01 


1.04±0.32 


0.11 


-0.03±0.02 


-0.40±0.61 


0.12 



Note. — Slopes, intercepts, and dispersions for each of the three samples at each X-ray energy, using 2500 A as the 
optical wavelength. Errors are at the la level. 
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Table 3. - ^bh Partial Correlation Coefficients 



X-ray Energy (keV) 


1500 A 


2500 A 


5000 A 


1.0 


-0.070 (1.84) 


-0.064 (1.68) 


-0.065 (1.71) 


1.5 


-0.036 (0.97) 


-0.030 (0.81) 


-0.032 (0.86) 


2.0 


-0.013 (0.35) 


-0.0071 (0.20) 


-0.0093 (0.26) 


4.0 


0.048 (1.37) 


0.053 (1.51) 


0.047 (1.38) 


7.0 


0.093 (2.66) 


0.096 (2.82) 


0.086 (2.53) 


10.0 


0.11 (3.14) 


0.10 (2.94) 


0.10 (2.94) 



Note. — The partial correlation coefficient listed is Kendall's ri2,3, 
with its significance (cr) in parantheses. Correlations with significance 
above 3.3cr (P > 99.9%) are considered to be significant. These are the 
results for the SPECTRA sample. 



Table 4. agx - ^hoi/^EM Partial Correlation Coefficients 



X-ray Energy (keV) 


1500 A 


2500 A 


5000 A 


1.0 


-0.13 (3.25) 


-0.12 (3.00) 


-0.12 (3.00) 


1.5 


-0.17 (4.25) 


-0.16 (4.00) 


-0.16 (4.00) 


2.0 


-0.19 (4.87) 


-0.19 (4.87) 


-0.19 (4.87) 


4.0 


-0.27 (7.11) 


-0.27 (7.11) 


-0.27 (7.11) 


7.0 


-0.32 (8.21) 


-0.32 (8.21) 


-0.32 (8.21) 


10.0 


-0.33 (8.25) 


-0.33 (8.25) 


-0.33 (8.25) 



Note. — The partial correlation coeflacient listed is Kendall's 

T12.3, with its significance (cr) in parantheses. Correlations with 
significance above 3.3cr (P > 99.9%) are considered to be signifi- 
cant. These are the results for the SPECTRA sample. 
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Table 5. a„T. - z Partial Correlation Coefficients 



X-ray Energy (keV) 


1500 A 


2500 A 


5000 A 


1.0 


-0.16 (4.32) 


-0.15 (4.05) 


-0.14 (0.036) 


1.5 


-0.11 (3.14) 


-0.097 (2.77) 


-0.087 (2.56) 


2.0 


-0.090 (2.65) 


-0.074 (2.18) 


-0.064 (1.88) 


4.0 


-0.040 (1.21) 


-0.027 (0.84) 


-0.021 (0.66) 


7.0 


-0.022 (0.67) 


-0.011 (0.33) 


-0.010 (0.31) 


10.0 


-0.015 (0.45) 


-0.013 (0.41) 


0.0086 (0.26) 



Note. — The partial correlation coefficient listed is Kendall's ri2,3, 
with its significance {a) in paranthcscs. Correlations with significance 
above 3.3cr (P > 99.9%) are considered to be significant. These are 
the results for the SPECTRA sample. 




Fig. 1. — The i band magnitude vs. redshift for three samples: CENSORED (circles), 
SPECTRA (triangles), and HIGHSNR (stars). 
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Fig. 2. — The aox — ^opt relation for the CENSORED sample. Blue circles indicate sources 
with S/N >6 (with a power-law fit to the X-ray spectrum), red circles indicate detected 
sources with S/N < 6 (no spectral fit, only X-ray fluxes), and black arrows indicate un- 
detected sources with X-ray upper-limits. The solid line is a linear leas t-squares fit to the 
data, with dispersions plotted as dotted lines. The dashed line shows the (jSteffen et al.ll2006l ) 
best-fit line. 
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Fig. 3. — (a)-(c) The — hpt relation for the SPECTRA sample, where the X-ray energy is 
1 keV and the optical luminosities are taken at 1500, 2500 and 5000 A. (d)-(f) The aox-hpt 
relation, where the X-ray energy is 10 keV and the optical luminosities are taken at 1500, 
2500 and 5000 A. Only sources with S/N > 6 are plotted. The colors signify redshift bins, 
which were used in determining the optical luminosities. Purple (z < 0.25), blue (0.25 < z 
< 0.43), light blue (0.43 < z < 0.55), green (0.55 < z < 1.8), orange (1.8 < z < 2.3), red 
(2.3 < z < 3.75) and black (z > 3.75). BAL and RL quasars are excluded, as are quasars 
with bad X-ray fits and quasars with X-ray fits that require significant intrinsic absorption. 
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Fig. 4. — The slope (a), intercept (b) and dispersion (c) of the aox — hpt relation with 
respect to the X-ray energy for the SPECTRA sample. The optical wavelength is held 
constant at 2500 A. The filled circles show the measured values, and the solid line shows 
the calculated values. Values are calculated by normalizing at the slope/intercept /dispersion 
for the relation at 2500 A and 2 keV, and then correcting for the baseline over which aox 
is defined. Measured values that lie away from the calculated baseline effect are affected by 
measurement errors and/or an intrinsic change in the — hpt relation. 
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Fig. 5. — The slope (a) and dispersion (b) of the aox — lopt relation for the CENSORED 
sample, as in Figure 3. 
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Fig. 6.— Optical color /^{g-i) vs. X-ray S/N for the CENSORED sample. The SPECTRA 
sample contains sources with S/N > 6, marked with a solid vertical line. X-ray weak sources 
(aox < -1.8) are marked as solid circles. 
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Fig. 7. — The slope (a) and dispersion (b) of the aox — Lpt relation for the HIGHSNR sample, 
as in Figure 3. 



0.4 



0.2 



^ 0.0 - 



CD 



o 








O 







o 








o 


^8 


o 


)CO 




o 



o o 




o 



o 



0.2 



■0.4 - 



-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 

log Lboi/LEdd 

Fig. 8. — Aaox vs. Lf,oi/LEdd- Defined at 4 keV and 2500 A, Aaox is the difference between 
the calculated aox and that predicted by the relevant aox — lopt relation. Colors show redshift 
ranges, as in Fig. 3. 
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Fig. 9. — Two notional quasar SEDs are shown, where specific luminosities are measured 
at 2500 A, 1 keV and 10 keV. Model (A) illustrates a source with lower lopt (10^^ ergs s~^) 
and (B) a source with higher lopt (10^^ ergs s~^). The X-ray luminosities {likeVi hokev) are 
determined from the appropriate — lopt relations, and the resulting imaginary power-laws 
defined by aox are shown as dashed lines. The average X-ray photon i ndex (< F > = 1.9) 
and the associated dispersion from the F— L^oi/L^Edd relation (cr = 0.33, iRisaliti et al.ll2009l ) 
are plotted as solid lines. 
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Fig. 10. — X-ray bolometric correction vs. ccox- The solid line is the best-fit, with the 
dispersion a{logK,) — 0.12 plotted as dotted fines. Only detected (S/N > 2), non-BAL, RQ 
sources with good X-ray fits are included. Where X-ray S/N > 6, sources with significant 
intrinsic absorption are excluded. 
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Fig. 11. — X-ray bolometric correction vs. Lhoi/LiEdd- The solid line is the best-fit, with the 
dispersion a{logn) = 0.38 plotted as dotted hnes. Sources included as in Figure [TOl 



